Abstract. We present the first numerical radiative trans- 
fer simulation of multiple light scattering in dust con- 
figurations containing aligned non-spherical (spheroidal) 
dust grains. Such models are especially important if one 
wants to explain the circular polarization of light, ob- 
served in a variety of astronomical objects. The optical 
properties of the spheroidal grains are calculated using the 
method of separation of variables developed by Voshchin- 



nikov & Farafonov (1993). The radiative transfer prob- 



lem is solved on the basis of the Monte Carlo method. 
, Test simulations, confirming the correct numerical imple- 
' mentation of the scattering mechanism, are presented. As 

■ a first application, we investigate the linear and circular 
polarization of light coming from a spherical circumstel- 
lar shell. This shell contains perfectly aligned prolate or 
oblate spheroidal grains. We investigate the dependence of 

, the results on the grain parameters (equivolume radius, as- 
(T) ' pect ratio) and the shell parameters (inner/outer radius, 

■ optical thickness). The most remarkable features of the 
~^ ' simulated linear polarization maps are so-called polariza- 

, tion null points where the reversal of polarization occurs. 
' They appear in the case when the grain alignment axis is 
perpendicular to the line of sight. The position of these 
t-H ' points may be used for the estimation of grain shape and 
, geometrical structure of the shell. The origin of null points 
' lies in the physics of light scattering by non-spherical par- 
, tides and is not related to the cancellation of polarization 
' as was discussed in previous models. The maps of circular 
polarization have a sector-like structure with maxima at 
O ' the ends of lines inclined to the grain alignment axis by 
£j ±45°. 
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1. Introduction 

It is now well established that the polarization of optical 
and near-infrared radiation from young and evolved stel- 
lar objects, reflection nebulae, and active galactic nuclei is 
mainly caused by dust grains. In many cases the observed 
polarization can be satisfactorily interpreted by light scat- 
tering on spherical grains. In particular, the orientation of 
the polarization vectors is often used for the investigation 
of the dust distribution and for identifying the location of 
the embedded illuminating source(s). 

However, even in the case of simple objects like ordi- 
nary reflection nebulae, deviations of polarization vectors 
from the direction perpendicular to a star were discov- 
ered more than three decades ago (Elvius & Hall 1967 ). 
These deviations are displayed in the outer filamentary 
parts of the Merope nebula and may be easily explained 
by light scattering on non-spherical grains aligned by a 
magnetic field along the filaments. Non-centrosymmetric 
polarization patterns have been observed in bipolar and 

19891 ), young stellar ob- 



cometary nebulae (Scarrott et al. 



jects (Hajjar & Bastien 1996| ), evolved stars (Kastner & 
Weintraub 1996| ), and are also clearly seen in polarization 
maps of comets (Dollfus & Suchail 1987). They may be at- 
tributed to non-spherical grains, although they could be 
caused by multiple scattering on spherical particles as well. 
Another effect which may be related to the light scattering 
by non-spherical grains is the wavelength dependence of 
the positional angle of polarization observed in red giants, 
AGB stars, a nd bi polar reflection nebulae (see, e.g., John- 
son & Jones 1991 ). The variations from blue to red may 
reach 20° — 60° and it is very difficult or even impossible 
to interpret this behaviour using spherical grains only. 

Recently, very high degrees of circular polarization of 
scattered light in the Orion molecular cloud were mea- 
sured by Chrysostomou et al. ( 2000 ). The authors sug- 
gest that the circular polarization is produced by aligned 



non-spherical grains. Multiple scattering of radiation by 
spherical or randomly oriented non-spherical grains re- 
sults in a much smaller circular polarization degree than 
observed. In addition, the interstellar polarization and 
polarized thermal emission phenomena prove that non- 
spherical grains exist in the interstellar medium. These 
effects arise because of dichroic extinction/emission of ra- 
diation by aligned non-spherical grains and were modeled 
with spheroidal grains (Voshchinnikov |l990j ; Kim & Mar- 



tin 



1995 Onaka 2000). 



Up to now, spheroidal grains have been used for the 
interpretation of scattered radiation only in the case of 
singl e sca ttering (Voshchinnikov |1998 ; Gledhill & Mc- 
Call 200C). The numerical simulation of the polarized 



radiation transfer through a medium with non-spherical 
grains, including multiple scattering, is extremely difficult. 
In the case of spherical grains nearly all previous simula- 
tions are based on Monte Carlo simulations. This method 
was used for the interpretation of polarimetric observa- 
tions of various stars (see, e.g., Daniel 1980 Voshchin- 



nikov & Karjukin 1994 ) as well as the production of po- 
larimetric maps of different ex tende d objects (e .g., B astien 
& Menard [1988^ Fischer et al. |l994| ; Wolf et al. |l998[ ). The 
most recent, major achievement of the Monte Carlo tech- 
nique — in respect to the solution to the radiative transfer 
(RT) problem — is the development of RT codes that al- 
low to calculate the spectral energy distribution of multi- 
dimensional dust configurations self-consistently (Wolf et 
al. [1999| ; Wolf & Henning p000| ) . 



Send offprint requests to: N.V. Voshchinnikov, nvvCDastro. 
spbu.ru 



The main aim of this paper is to provide the first ra- 
diative transfer simulations including scattering, extinc- 
tion, absorption, and re-emission of radiation by aligned 
spheroidal dust grains. Our formalism will be applied to 
several simple model configurations. Because of the com- 
plexity of the problem, we decided to consider the separate 
mechanisms step by step. In this paper, the basic theory 
is presented and the main numerical features of the sim- 
ulation are described (Sect. §). Furthermore, we consider 
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multiple scattering of light by spheroidal grains and the 
resulting linear and circular polarization (Sect. [?]and|4|). 

2. Basics 

2.1. Dust grain parameters 

Let us suppose that the dust grains are prolate or oblate 
homogeneous spheroids with an aspect ratio a/b, where 
the quantities a and b are the major and minor semiaxes of 
a spheroid, respectively. The particle size is characterized 
by the radius ry of a sphere with the same volume as that 
of the spheroid (equivolume radius). The major semiaxis 
a of the spheroid is connected with ry by 

.-„.(!)«■ (i) 

for prolate spheroids and 

a = rv( ¥ ) (2) 
for oblate spheroids. The optical properties of a particle 
also depend on the complex refractive index of the grain 
material m\ and the angle a between the rotation axis of 
the spheroid and the wave-vector direction (direction of 
the incident radiation; 0° < a < 90°, see Fig. |l|). 

2.2. Radiative transfer concept 

In this section we outline the main ideas of our solution 
of the RT problem taking into account multiple scatter- 
ing and dichroic absorption, and reemission by spheroidal 
grains. Our solution is based on the Monte Carlo method 
which is described only briefly here. For a more general 
overview about RT based on the Monte Carlo method, 



we refer to the papers of Wolf ct al. (1999) and Wolf & 
Henning (2000) (see also Fischer et al. 1994). 

In our solution to the RT problem, the radiation en- 
ergy is partitioned into so-called weighted photons (in 
the following we will use the term "photon"). The en- 
ergy, intensity, and polarization of a weighted photon 
are described by its wavelength A and the Stokes vector 
/ = (I, Q, U, V) T . The transformation of the Stokes vector 
by the ith scattering can be described by a 4 x 4 (Muller) 
matrix F: 
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The main features distinguishing the light scattering 
by aligned non-spherical particles and spherical particles 
are: 

1. Azimuthal dependence of the scattered radiation; 

2. Linear polarization in the forward and backward di- 
rections; 

3. Deviation of the positional angle of linear polarization 
after first scattering from the direction perpendicular 
to the illuminating source and 

4. Circular polarization after first scattering. 



These features result from the fact that in the general case 
of non-spherical particles — in contrast to spherical grains 
- all 16 elements of the scattering matrix are non-zero 
(see Voshchinnikov & Farafonov 1993). 

The spatial distribution of the illuminating star(s) and 
the dust grains is defined inside a convex model space. The 
free path length I of a photon (distance between the point 
of emission and the first point of scattering/absorption or 
two points of scattering/absorption) is determined by the 
dust density distribution nd(Ri) at the coordinate Ri and 
the extinction cross-section for the transmitted radiation 

Ccxt(Ri) 



i=l 



n d (Ri) ■ C oxt (Ri) ■ Ah 



(4) 



The quantity i cn( j represents the number of integration 
steps along the free path length, and the quantity AZj is 
the corresponding geometrical step width. The free mean 
path length 



i = J2 Ah 



(5) 



can be derived from the optical depth f ex t as follows 
f oxt = -ln(l-C) ■ ' (6) 

Here, £ is a random number uniformly distributed on the 
interval [0, 1] (see Wolf et al. |l999| for details). 

The extinction cross-section C ox t is calculated taking 
into account the state o f pola rization of the incident radi- 
ation (see, e.g., Martin 1974 ): 



<--cxt — Xl^cxt 
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Here, $ is the angle between the particle frame and the 
laboratory frame. (Io,Qo,Uq,Vo) is the Stokes vector of 
the incident light. The extinction cross-sections C, 



TE.TM 



are connected with the efficiency factors Q 



TE.TM 



TE.TM 



GQ 



TE.TM 



by 



(8) 



where the quantity G is the geometrical cross-section of 
the dust grain (the area of the particle shadow) . The su- 
perscripts TM and TE are related to the two cases of 
the polarization of the incident radiation (TM and TE 
modes). The geometrical cross-section G can be derived 
as follows 



G(a) 



U) • w sm a+ 



' V 



cos 2 a 



1/2 



(9) 



cos 2 a 



-l 1/2 



(10) 



TM,TE 



for a prolate spheroid and 

f a \ 2/3 / a \ — 2 
G(a) = 7rr v • (^-J • ^-J sin 2 a + 

for an oblate spheroid. The efficiency factors Qoxt ' ~ are 
calculated from the solution to the light scattering prob- 
lem for spheroids obtained by the method of separation 
of variables (see Voshchinnikov & Farafonov 1993 for de- 
tails). 

We would like to note that the extinction cross-section 
Cext (see Eq. (Q)) has to be determined at each integration 
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Scattered 
ray 



Incident ray 

Fig. 1. Scattering geometry for prolate homogeneous 
spheroid. The origin of the coordinate system is at the 
centre of the spheroid while the z-axis coincides with its 
axis of revolution. The angle of incidence a is the angle in 
the x — z-plane between the direction of incidence and the 
z-axis. The scattered field in the far-field zone is repre- 
sented in a spherical coordinate system by the angles 9, if 
in the reference system related to the spheroid's rotation 
axis (particle frame) or by the angles O, $ in the reference 
system related to the scattering plane (laboratory frame). 

step along the photon path (see Eq. (Q)). In the case of 
the RT in a medium with aligned non-spherical grains, the 
step width along the photon free path is determined not 
only by the scale of the density structure (as for spherical 
particles) but also by the degree of spatial variations of 
the grain alignment. Owing to this and the more complex 
treatment of the scattering process compared with spher- 
ical grains, the total computation run-time is by an order 
of magnitude larger for aligned spheroids than for spheres. 

Finally, depending on the model geometry and the 
dust density distribution, the photon leaves the model 
space after a certain number of interactions with dust 
grains (scattering, absorption). Then, it will be "observed" 
by array-detector-like planes arranged around the model 
space. The Stokes vector has to be projected onto the re- 
spective plane. For a detailed description of the projection 
formalism we refer to Fischer et al. ( |1994| ). 

2.3. Scattering 

The solution to the problem of light scattering by 
spheroidal particles is given in a particle reference sys- 
tem related to the rotation axis of a spheroid (in spherical 
coordinates 9, ip), while the RT is considered in the labo- 
ratory reference system (in spherical coordinates O, $; see 
Fig. [y). The probability of a photon of being scattered into 




Fig. 2. The linear polarization of radiation P\(Q) after the 
first scattering by a prolate spheroidal particle with refrac- 
tive index m = 1.7 + 0.03i (typical of astronomical silicate 
at optical wavelengths), equivolume radius ry = 0.20 /xm, 
aspect ratio a/b = 4 and the incidence angle a — 0°. The 
wavelength of incident radiation is A = 0.628 /im. To avoid 
skipping of local minima/maxima a step size smaller than 
5° for is required. 



the direction (0, $) is proportional to the value of the first 
element of the scattering matrix Fn(a,Q,&). We trans- 
formed the scattering matrix from the particle to the lab- 
oratory frame. Then, the scattering matrix is calculated 
using the standard expressions (Bohren & Huffman 1983). 
In order to check the transition from the particle system to 
the laboratory one, the transformation was also done for 
light scattering by spheroidal particles in the quasistatic 
approximation where the solution can be written in both 
coordinate systems (see Farafonov et al. 



2001 



for details). 

The elements of the scattering matrix were normalized 
using the standard procedure (Bohren & Huffman |1983 ) 

——?—-—[ [ F n («,©,$) sin 0d0d$= 1(11) 

K l°sca ' °sca > JO JO 

where C sca is the scattering cross-section and k = 2n/X 
the wavenumber. The correctness of these computations 
was verified using the relationships for matrices, describ- 
ing scattering by small particles (see Hovenier & van der 
Mcc 200C and references therein). 

The calculation of the elements of the scatter- 
ing matrix F(m, ry, a/b, a, 0, $) and efficiency factors 
Qscn cxt ( m i r v, a/b, a) for a single particle with the fixed 
equivolume radius ry , the aspect ratio a/b, and the refrac- 
tive index m requires several minutes of computer timcQ. 
These quantities were precalculated before starting the RT 
modelling. In order to guarantee that the main features 
(local minima/maxima) of the scattering matrix elements 
(and/or derived quantities, for instance, the polarization 
of singly scattered light) are present in the discrete dis- 
tributions as a function of a, 0, and $ (see Figs.|| and 



A Pentium II 450 MHz processor was used for our calcula- 
tions. 
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Fig. 3. The degree of linear polarization in the midplane 
of the polarization map shown in Fig.|j[A] in dependence 
of the projected distance. Three different discretizations 
for the angles a, O, and <& were applied (Aa = AG = 
A$ = 2°, 10°, and 15°). Only the smallest step width 
gives a sufficient approximation to the analytical solution 
presented in Fig.|| [Bl] (see Sect. 2.2 for discussion). 



Before applying the albedo matrix A part , it has to be 
rotated into the particle frame using the rotation matrix 
R 

(\ 0^ 

cos 2* sin 2* 



R = 



-sin 2* cos 2* 







(16) 



1/ 



The absorption process can then be described by the ex- 
pression 

I after = R Ap ar t. -R /before , (17) 

whereby /before (/after) is the Stokes vector of the radiation 
before (after) absorption. 

In the case of (re)emission of radiation by a spheroidal 
particle, the Stokes vector of the photon can be written 
as follows (particle frame): 
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B x (T d ). (18) 



part. 



Here, C a b s is the absorption cross-section, B\(T<i) the 
Planck function and Td the grain temperature. 



3. Test simulations 



||), we use step widths of Aa = AG = A$ = 2° in our 
simulations (see Fig. ||). 



2-4- Absorption and Reemission 

The absorption is described by the albedo A of the par- 
ticle which is defined as the ratio of the scattering and 
extinction cross sections. It is equal to 



A 



TM,TE 



(to, ry, a/b, a) 



a 



TM.TE 



(m, ry, a/b, a) 



c cxt {m,r v ,a/b,a) 



(12) 



(13) 



for totally polarized light and 

r«TM _|_ /-VTE 

AW{m,rv,a/b,a)= 

^ext "T" u ext 

for non-polarized light. 

While in the case of absorption by spherical particles, 
the Stokes vector / can be multiplied by a scalar value 
of the albedo, for spheroidal particles the following 4x4 
matrix is needed (in the particle frame): 

/ h h \ 

h h o 



A 



part. 
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o o h o 
V o o o h J 



(14) 



part. 



Here, the quantities l\ and I2 are defined as 
h = A 

h 
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(15) 



The best way to test any RT code is the compari- 
son of its results with those obtained with another RT 
code/technique or an analytical solution. Because the RT 
code presented here is the first that provides the treat- 
ment of multiple scattering of polarized light by aligned 
non-spherical particles, a lack of suitable test cases for 
comparison is evident. Nevertheless, in the following we 
present test calculations which confirm the correct imple- 
mentation of the new scattering mechanism into our RT 
code. 

For test simulations we used particles with the follow- 
ing parameters: 

— shape of the grains: oblate, 

— refractive index: m = 1.7 + 0.03i, 

— equivolume radius: ry — 0.1 /im, 

— aspect ratio: a/b = 2.0, 

— wavelength of incident radiation: A = 0.628 /im. 
Then the grain size parameter is xy = 2itry/\ = 1. 

3.1. Single scattering in a geometrically and optically 
thin hemisphere 

Firstly, we consider single scattering in a thin layer of 
oblate spheroids at the surface of an hemisphere (see 
Fig. |). The ratio of the inner radius (n n ) of the shell 
to the outer one (r ou t) amounts to nn/fout = 0.99. The 
non-rotating particles show a static picket fence orienta- 
tion with their rotation axes being aligned parallel to the 
z-axis of the global coordinate system of the model space. 
The density profile n(r) in the shell follows a power law: 
n(r) oc r^ 1 . (19) 
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Fig. 4. Model configuration for comparison with the re- 
sults of Matsumura & Seki ( 1996 ). It consists of a thin 
dust layer at the surface of a hemisphere. The ratio of the 
inner radius to the outer one amounts to ri n /r ou t = 0.99. 
The shell contains oblate grains, whereby their rotation 
axes (b) are oriented parallel to the z-axis. The dark grey 
oval pattern symbolizes the aligned grains. The configu- 
ration is shown for y = (irz-plane, the azimuthal angle 
$ = 180°). An observer sees the radiation at the scattering 
angles 9 = 0°-90°. In this case the condition Q + a = 90° 
is satisfied in all points of the dust layer. 



The optical thickness of the shell at the wavelength of 
the simulated RT (A = 0.628 fim) is 0.2 along the z-axis, 
i.e. r z (0.628/im) = 0.2. Because of the chosen grain align- 
ment, r z is proportional to the extinction cross-section in 
the case of light incident parallel to the grain rotation axis, 
tq = t z rx C cx t(0°). In other directions, the optical thick- 
ness of the shell may be found as r a = T z C ex t{a)/C ex t(0°) 
and, in particular, along the line of sight (x-axis) and the 



y-axis we have r x = r y = r gQ o = r z C oxt (90° )/C ex t(0°). Us- 
ing Eq. [?] for non-polarized stellar radiation and the stan- 
dard model described above, we have r go o w 0.35 (prolate 
grains) and r gQ o w 0.13 (oblate grains), if r (0.628 fim) = 
0.2. For more elongated or flattened particles (a/6 = 10), 
these values become equal to t qq o w 0.68 and 0.082 for 
prolate and oblate grains, respectively. 

In contrast to polarization maps obtained for model 
configurations with the same geometry, but containing 
spherical dust grains, the polarization pattern is not 
centro-symmetric (see Fig. ||[A]). This is because of the 
azimuthal asymmetry of light scattering by aligned non- 
spherical particles. The net polarization in this case is not 
equal to zero. We compare our results with those obtained 
by Matsumura & Seki fll996p pl In Fig. | the comparison 
of the polarization degree as a function of the projected 
distance from the centre along three different lines in the 



polarization map is shown. These graphs as well as the 
polarization pattern (see for comparison Fig. 2 from Mat- 
sumura & Seki 1996) are in perfect agreement. 



Qualitatively, the polarization pattern is quite similar 
to that of a young stellar object surrounded by a spher- 
ical shell and an optically thick disk seen edge-on. There 
the parallel polarization pattern is due to light scattering 
on the disk surface while the centro-symmetric polariza- 
tion pattern is due to light scattering in the circumstellar 
environment. Corresponding RT simulations with spher- 



ical dust grains hav e bee n performed by Menard (1988) 
and Fischer et al. ( 1994 ). The most interesting feature 



2 The nume rical data in tabular form were provided by Mat- 
sumura (199£ ). 



of the pattern is the presence of polarization null points 
(see Figs. |][A] and |][B1]). They have been obser ved in a 
number of objects (see discussion in Sect. 4.2.3 ). Some- 
times, these null points are interpreted as the indicators 
for the "outer edge" of a circumstellar disk (transition 
from the opticall y thi ck disk to the optically thin shell; 
see Fischer et al. 1994 ). Although this explanation is very 
plausible, the simulations presented here show that uni- 
formly aligned grains (possibly by an external magnetic 
field) may produce a similar polarization pattern. A de- 
tailed discussion of the polarization null points as a func- 
tion of the model parameters is given in Sect. ||. 

3.2. Multiple scattering 

To test the RT with multiple scattering, we consider model 
configurations for which the dust density distribution, the 
dust grain orientation, and the irradiation characteristic 
are centro-symmetric with respect to the observer's posi- 
tion. The resulting polarization pattern is therefore also 
centro-symmetric: the polarization vectors are oriented 
perpendicular or parallel to the radius vector of the con- 
figuration, in dependence of the refractive index and the 
grain size parameter xy. Thus, the integral (net) linear 
and circular polarization have to be zero. Because this 
property of centro-symmetric configurations does not de- 
pend on its optical thickness, it can be used as a test of 
the single and multiple scattering case. 

We considered the following dust configurations con- 
taining cither oblate or prolate grains: 

1. Spherical dust shell with a point-like radiation source 
in the centre and particle rotation axes either parallel 
or perpendicular to the radius vector; 

2. Cylindrical dust shell with a point-like radiation source 
on its rotation axis and particle rotation axes cither 
parallel or perpendicular to the radius vector which 
is oriented perpendicular to the rotation axis of the 
cylinder. 

While for the first configuration the net polarization 
is equal to zero for every position of the observer, for the 
second configuration this is the case only if the observer's 
line of sight is oriented perpendicular to one of the two 
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Fig. 5. The results of test calculations for the comparison with the data of Matsumura & Seki ( 1996| ). [A]: Polarization 
map of the configuration shown in Fig. |]. Only single scattering is considered. [B1]...[B3]: The degree of linear 
polarization as a function of the projected distance from the centre along three different directions. The projected 
distance r/r out is connected with the scattering angle as r/r out — arcsinO. [Bl]: Pi in the midplane of polarization 
map; [B2]: Pi perpendicular to the midplane; [B3]: P along a line being inclined by 45° to the midplane. 



cylinder footprints. In all cases we arrived at the expected 
results. 



4. Results and discussion 



In Sect. 4.1 j light scattering in a geometrically thin hemi- 



sphere and sphere and a geometrically thick spherical 
shell is described. This preliminary modelling provides 
a "smooth" transition from the test models discussed in 
Sect. H to more realistic dust configurations. The basic 



modelling is presented in Sect. 'L2, where a geometrically 
thick shell is considered and the resulting polarization as 
a function of the main model parameters is investigated. 

If one considers the RT in a medium consisting of 
spherical dust grains, scattering is the only mechanism 
to produce or to modify the polarization of the incident 
radiation. In the case of aligned spheroidal grains, the 
processes of photon absorption and re-emission may also 
fabricate or change polarization. Because of the obvious 
complexity of the problem, we separated the consideration 
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Fig. 6. The degree of linear polarization in dependence 
of the projected distance for tilted dust configurations. 

of these three mechanisms of the polarization formation. 
The aim of the simulations discussed below was to outline 
the basic features of polarization arising from the light 
scattering by spheroidal grains. The additional influence 
of absorption and re-emission will be discussed in future 
publications. 

For this reason, we apply the following approxima- 
tions: 

1. The grain albedo is assumed to be that of spherical 
grains with the equivolume radius ry and is therefore 
a scalar. Then, the absorption process decreases the 
intensity of the incident radiation but does not modify 
the state of polarization. 

2. We do not consider thermal re-emission of radiation. 

If not stated differently, for the simulation presented 
in this section the dust grain and model parameters as de- 
scribed in Sect. |||have been used. The observer's position 
is inclined by i — 90° to the z-axis (r z (0.628 /an) = 0.2) 
of the global coordinate system of the model space (see 
Fig. |). 

4-1- Simple models 

Let us consider again the geometrically thin hemispherical 
model. While the polarization pattern in Fig. || is shown 
only for a position perpendicular to the z-axis (i = 90°), 
Fig. H demonstrates how the behaviour of the degree of lin- 
ear polarization in the midplane of the polarization map 
appears for different inclinations of the z-axis to the line 
of sight. It turns out that the local polarization minimum 
(the "polarization null point" or the reversal of polariza- 
tion in the case i = 90°) are present at large inclination 
angles i only. It is important to note that the polariza- 
tion pattern (and therefore also the net polarization de- 
gree) for aligned non-spherical particles depends on the 
viewing angle even in the case of a one-dimensional den- 
sity distribution. In contrast to this, no dependence of the 




r y /r out 



Fig. 7. The degree of linear polarization in dependence 
of the projected distance. Influence of the backward scat- 
tering is illustrated. 



polarization on the inclination angle exists for spherical 
particles. 

As a first step away from the very artificial configura- 
tion described above, we consider the RT in a thin spheri- 
cal shell instead of the hemisphere shown in Fig. ||. Then, 
backward scattering in the hemisphere opposite to the ob- 
server's position (where = 90° + a) has to be taken into 
account. The results are shown in Fig. 0. The addition 
of the backscattered radiation shifts the position of the 
minimum of the linear polarization to a smaller projected 
radius. 

In Fig. ||, the circular polarization for this model of a 
geometrically and optically thin dust shell is shown. The 
circular polarization is equal to zero in the midplane and 
along the line through the centre perpendicular to the 
midplane. The net circular polarization from the config- 
uration is also equal to zero. A counterpart of the lin- 
ear polarization null points in the midplane could not be 
found. The maximum circular polarization is observed at 
the ends of lines inclined to the z-axis by ±45°. These 
points correspond to the light scattering by particles with 
values of a = 45°, 9 = 90° and $ = 90° (270°). 

It turned out that the pattern of the linear polariza- 
tion strongly depends on the inner radius of the spherical 
shell. This is demonstrated in Fig. [| where the degree 
of linear polarization in the midplane is shown for three 
ratios of the inner to the outer radius of the shell. Decreas- 
ing this ratio (and therefore decreasing the relative size of 
the inner cavity around the illuminating star), the char- 
acteristic polarization null points shift towards the star. 
It has to be noted, that generally the polarization null 
point is not located at the position of the projected in- 
ner radius. The location, relative depth and width of the 
area where the polarization changes the sign depends on 
the geometric al thic kness of the shell and grain parame- 
ters (see Sect. 4.2.2 ). The power- law density profile used 
(see Eq. (|l9|)) leads to a preferred photon scattering near a 
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Fig. 8. [A] The circular polarization of the scattered radiation arising from the configuration shown in Fig. [l| 
Deviations from the mirror symmetry result from the statistical noise of the Monte Carlo simulation. [B] The same 
as grey plot to emphasize the symmetry of the spatial dependency of the circular polarization degree. 
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Fig. 9. The degree of linear polarization in dependence of 
the projected distance. The effect of the transition from a 
geometrically thin to a geometrically thick spherical shell 
is illustrated. 



narrow region at the inner boundary of the shell. The shift 
of the polarization null point is therefore mainly caused 
by the change of the inner radius ri n . 

4-2. Envelope models 

Now we will consider geometrically thick spherical shells. 
To represent such models, the ratio of the inner radius of 
the shell to the outer radius is fixed at r- m /r ou t = 0.5. 



4.2.1. Transition from optically thin to optically thick 
shell 

In the case of a small optical thickness, the main contri- 
bution to the resulting polarization is from first or sec- 



ond scattering events (Voshchinnikov & Karjukin 1994 



Voshchinnikov et al. 1995). As one can see from Fig. h0 



the polarization P\(r y /r ou t) in the midplane decreases 
with increasing optical thickness. The same behaviour was 
found for the polarization perpendicular to the midplane 
P\{r z /r on t). In general, the position of the local polariza- 
tion minimum is not strongly influenced by the optical 
thickness. 

Integrating the images for the Stokes parameters 
/, Q, U, we estimate the total linear polarization from our 
spherical dust configuration as one observes from the posi- 
tion perpendicular to the z-axis. Because of the symmetry, 
it is the same along the x and y-axes 



sea. 90 



1,90" 



F*exp 



90 



F 



(20) 



where F 



sea, 90 

and F p are the scattered and polar- 
sea, 90 scaj go 

ized scattered radiation from the configuration. For the 
models presented in Fig. |l0|, the resulting polarization is 
Pi = 4%, 9%, and 8%, if r z = 0.2, 1 and 2, correspondingly. 
Excluding the stellar radiation (first term in the denom- 
inator in Eq. (pfj|)), we can find the maximum degree of 
the observed linear polarization 
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Fig. 10. The degree of linear polarization in dependence 
of the projected distance. The effect of the transition from 
an optically thin to an optically thick spherical shell is 
illustrated. Note that at this Figure the resolution along 
the y-axis is much smaller than in comparison with other 
Figures. 



It is equal to P. 



1,90 



25±1%, 15±1%, and 9.1 ±1%, 



if r z = 0.2, 1 and 2, correspondingly. 

4.2.2. Dependence on dust grain parameters 

Figure O shows the midplane polarization Pi(r y /r out ) for 
oblate and prolate grains of three different equivolume 
radii (ry = 0.05 /im, 0.1 /im, 0.2/im). The behaviour of 
polarization for particles with ry = 0.05 /im and 0.1 /im 
is rather similar, while the polarization produced by the 
largest particles (ry = 0.2 /im) is the smallest and is al- 
most independent of the projected radius. This feature 
reflects the fact of decreasing polarizing efficiency of large 
particles (Voshchinnikov et al. 2000] ). Because the man- 
ner of polarization for prolate grains in the direction per- 
pendicular to the midplane is similar to that for oblate 
grains in the midplane and vice versa (see also Fig. |l3| ), 
the resulting polarization for particles of both types is 
quite similar. The net polarization — excluding the direct, 
non-scattered stellar radiation — amounts to P i gQ o max = 
20% (16%), 25% (20%) and 11% (10%) for oblate (pro- 
late) grains with equivolume radii ry = 0.05 /im, 0.1 /im 
and 0.2 /im, respectively. 

The distribution of the linear polarization P\(r y /r ou t) 
in the midplane of the spherical shell is plotted in Fig. [l2] 
for particles with four different aspect ratios a/b (a/b = 
1.1, 2, 4, 10). The polarization minima in the midplane 
(r z = 0) appear for oblate grains only. Their position 
shifts to outer projected radii r y /r ou t with increasing as- 
pect ratio. If a/b — ► 1, the location of the polarization null 
point leads to r y /r ont — * 0. This is consistent with the be- 
haviour of the degree of linear polarization for spherical 



grains where the null polarization can be found at the 
centre for the scattering angle 9 = 0°. The resulting po- 
larization for oblate and prolate grains of different shape 
is given in Table [l]. Note that the polarization produced 
by oblate particles is larger in all cases. 



Table 1. The net linear polarization from the spherical 
shell containing perfectly aligned spheroidal grains (m = 
1.7 + 0.03i, ry = 0.1 /im) as a function of the grains' axes 
ratio a/b. The direct stellar radiation is excluded. 



a/b 


Oblate particles 


Prolate particles 


1.1 


3.6% 


3.0% 


2 


25% 


20% 


4 


43% 


31% 


10 


57% 


41% 



The intensity maps of our configuration (i — 90°) show 
an additional effect. In Fig. [l3|, the images with overlaid 
polarization pattern, corresponding to the intensity and 
linear polarization of scattered light, are shown for oblate 
and prolate grains with the aspect ratios a/6=l.l and 2. 
In the case of nearly spherical grains (a/b — 1.1), the ring- 
like structure is clearly seen for oblate as well as for prolate 
grains. Otherwise, in the case of a/b = 2 the intensity max- 
ima appear at the points of intersection with the z(y)-axis 
for oblate (prolate) grains. In both cases, the correspond- 
ing axis is parallel to the minor (b) axis of grains. Increas- 
ing the aspect ratio a/b, this effect is strengthened. Note 
that the net polarization from the spherical configuration 
arises in the inner parts of the shell because the polariza- 
tion pattern from the outer regions is centro-symmetric 
and the polarization would be cancelled. The direction of 
the resulting net polarization is parallel to the major (a) 
axis of grains that is similar to the case of dichroic polar- 
ization. 

4.2.3. Null points on polarization maps as indicators of 
scattering by non-spherical grains 

The model discussed in this paper is still too simple in 
order to make a detailed comparison with observations. 
However, we will discuss qualitatively the behaviour of 
null points produced by scattering on non-spherical grains. 

Polarization null points were observed in a variety of 
objects ranging from reflection nebula and circums tellar 
structures to e xtern al galaxies (e.g., Scarrott et al 



1989 



Sc arrott et al. |l990j ; Gledhil l fc S carrott |l989j ; Asse lin et 
al. 1996 ; Chrysostomo u et a l. |l996| ; Hajjar fc Bast ion 1996 ; 
Kastn er fc W eintraub |l996| ; Wood & Jones |1997| ; Gledhill 
et al. 200l| ). They appear approximately symmetrically 
on either side of the central illuminating source^]. In all 

3 It should be mentioned that the measurable polarization 
— especially near the position of polarization null points — 
sensitively depends on the spatial resolution of the polariza- 



Wolf et al.: Multiple scattering of polarized radiation by non-spherical grains 

Oblate grains Prolate grains 

100 

If [%] 



11 




0.6 r y /r ou t 




0.6 r y /r ou t 



Fig. 11. The degree of linear polarization in dependence of the projected distance. The effect of the variations of the 
equivolume radius ry for oblate and prolate grains is illustrated (ry = 0.05 /^m, 0.1 /im, 0.2 ^m). 
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Fig. 12. The degree of linear polarization for the case of oblate and prolate grains in dependence of the projected 
distance. The effect of the variations of the grain shape (the aspect ratio) is illustrated (a/b = 1.1, 2, 4, 10). 



cases the origin of null points was attributed to the com- 
bined action of two polarization mechanisms, where the 
first (initial or internal) mechanism produced linearly po- 
larized radiation whereas the second (external) mechanism 
led to the compensation of this polarization. Three pos- 
sible models were qualitatively discussed by Scarrott et 
al. ( |1989| ): 

(i) Polarization due to scattering by spherical grains 
and its cancellation via dichroic extinction in an external 
medium consisting of aligned non-spherical particles; 

(ii) Polarization due to scattering by spherical grains 
in an optically thick internal medium and its cancellation 
after the next scattering by spherical grains in an external 
optically thin medium; 

(iii) Polarization of light from the central source due 
to an unspecified mechanism and its compensation due to 



the second scattering by spherical grains in an external 
medium. 

These models could be called "hybrid models" because 
they require two steps to produce the observed polariza- 
tion null points. The model (i) is rather naturally applied 
to the modelling of polarization maps of galaxies (Wood & 
Jones 1997). The multiple scattering model (ii) was used 
by Bastien & Menard ( |1988| ) and Fischer et al. ( |1994[ ) for 
interpretation of polarization maps of young stellar ob- 
jects. In this case, the centrosymmetric polarization pat- 
tern is obtained in outer optically thin parts of the shell, 
whereas the aligned polarization vectors appear in the 
dense inner parts of this shell. The polarized source model 



(iii) was discussed by Notni (1985) and Gledhill (1991) 



tion maps. Therefore, it is more likely to find the polarization 
minima but not real null points. 



Using this model, it is possible to explain the polarization 
pattern observed i n sev eral bipolar and cometary nebu- 
lae (e. g., G ledhill 1991 ), but as it was shown by Clark 
et al. (2000) that the models with spherical grains fail to 
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produce even the 5 % circular polarization found in the re- 
flection nebulosity around the Herbig Ae/Be star R CrA. 
In this case, the authors suggested scattering by aligned 
non-spherical grains as the operating mechanism. 

In contrast to the hybrid models, different polarimetric 
effects can arise as a result of light scattering by aligned 
non-spherical particles. In our simple model, the location 
of polarization null points depends on the particle shape, 
size and the inner to outer radius ratio (see Figs. |9[ |ll| and 
[l2] ). The wavelength dependence of the location of null 
points is not very strong and their position only slightly 
approaches the central (stellar) position when going from 
the V to K band as it follows from Fig. 14 A simi- 



frill 
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0.44 
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Fig. 14. Projected positions of null points on polarization 
maps in dependence of wavelength and their statistical 
errors for the case prolate grains with aspect ratio a/b = 2 
(the similar model is shown in Fig. right lower panel). 
The calculations were performed for effective wavelengths 
of Johnson photometric system. 



lar conclusion can be reached from the observed polar- 
ization maps of the v ariable Hubble's nebula NGC 2261 
(Scarrott et al. 1989 ) and the infra red reflection nebula 

The results of As- 



199G) 



GSS 30 (Chrysostomou et al. 
selin et al. (1996) for V376 Cas show that the position of 
the null points is also shifted as a function of the binned 
box size used in averaging the polarization. However, it 
should be noted that near the null points the observed 
signal is usually very small and the errors are large. 

Finally, we can indicate that the scattering by almost 
spherical grains is marked by a ring-like structure in the 
intensity maps whereas the null point are located inside 



4 In the case of polarization null points which have been ob- 
served at the location of the outer boundary of circumstellar 
disks around young stellar objects (transition from the opti- 
cally thick disk to the optically thin shell), the same behaviour 
is expected since this boundary is slightly shifted inward with 
increasing wavelength. 



this structure rather close to the stellar position (Fig. [L3L 
left side). In the case of aligned non-spherical particles 
(Fig. 0, right side) , the null points at the inner boundary 
of the spherical shell arc accompanied by symmetric arc- 
like structures in the intensity maps. There is some simi- 
larity of such maps with those of small reflection nebulae 
published by Gledhill et al. (2001): the "spherical case" 
is seen in the source IRAS 19114+0002 while the "non- 
spherical case" occurs in the maps of IRAS 22223+4327, 
IRAS 19500-1709 and IRAS 17106-3046. 

We should emphasize that in the presented model the 
origin of null points is caused by light scattering by non- 
spherical particles, but is not related to the cancellation 
of polarization as in hybrid models. 



5. Conclusions 

The paper contains the first solution to the RT problem 
including multiple light scattering in dust configurations 
with aligned non-spherical dust grains. Our main aim was 
to discuss the effects of light scattering by spheroidal dust 
grains in a simple (spherical) dust configuration, contain- 
ing perfectly aligned particles. 

The most remarkable features of the simulated linear 
polarization maps are the polarization null points where 
the reversal of polarization occurs. They appear when the 
grain alignment axis is perpendicular to the line of sight. 
Symmetrically to the polarization null points, we found 
maxima in the intensity maps. 

In contrast to spherical grains, even single scattering 
by spheroidal grains may cause circular polarization. The 
maps of circular polarization have a sector-like structure 
with polarization maxima at the ends of a line inclined to 
the grain alignment axis by ±45°. 

Based on our theoretical investigations of light scatter- 
ing, dichroic absorption, and (re)emission by spheroidal 
grains presented in Sect. ^[ the next steps of the numeri- 
cal radiative transfer simulations will include the consid- 
eration of light scattering by partly aligned rotating par- 
ticles and the calculation of the polarized re-emission of 
grains in different dust configurations. This will imme- 
diately lead to an improved interpretation and therefore 
better understanding of mid-infrared and submillimeter 
polarization where re-emission and dichroic absorption by 
partly aligned non-spherical dust grains are the main po- 
larization mechanisms. Thus, it provides a basis for inves- 
tigations of the structure of the magnetic fields in these 
objects which are in most cases assumed to cause grain 
alignment. 
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Fig. 13. Intensity maps with overlaid polarization pattern of a spherical shell containing perfectly aligned dust grains 
(for the description of the model configuration see Sect. 4.2.2). Due to chosen shell and grain parameters, the brightest 
knots appear at the inner boundary of the shell. 



